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Abstract 

The coupling constant of p — * 777 decay is calculated in the framework of light 
cone QCD sum rules. A comparison of our prediction on the coupling constant with 
the result obtained from analysis of the experimental data is performed. 
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1 Introduction 



Radiative transition between vector(V) and pseudoscalar(P) mesons represent an important 
source of information on low energy hadron physics. These transitions are governed by 
magnetic dipole (Ml) radiation of the photon and had played one of the central roles 
for checking the predictions of quark model and SU(3) symmetry, as well as it was very 
useful in the determination of the magnetic dipole moment of iV*(1535) in 7^ — > r)N 
process [3]. Recently the theoretical activity on the VP7 magnetic dipole transition have 
increased(see [1] and references therein), in particular, due to the fact that the analysis of 
the radiative V — > P7 decay with 77 and rj mesons in final state can provide insights to the 
long standing issue of the 77 and rj mixing (for review see [2] and references therein). The 
usual parametrization of r\ — rj mixing in the octet-singlet basis, which will be used in this 
work, is as follows: the current-particle matrix elements are defined as 

<0| J" |P(p)> = = 8; 0; P = V , rf) (1) 

where Jf , the SU(3)f octet axial vector current, is given by 

J ln = ^("7^75" + d>y^ 5 d - 2s-f^ 5 s) (2) 
and J° M , the SU(3)f singlet current is given by: 



= -^( M 7^75« + d^ 5 d + S7^7 5 s) (3) 



Two mixing angles #8 and #0 are required in order to consistently describe mixing [4]. 
Accordingly the couplings in Eq.(l) can be defined as follows 

ft = f 8 cos9 8 , f° = -f sin9 (4) 
f*=f 8 sin0 8 , f°=f cos9 (5) 

Alternatively, two independent axial vector currents with distinct flavors can be consid- 
ered 

J 5» = S7 M 7 5 s (6) 
and the couplings of these currents with 77 and 7/ mesons are defined similarly to Eq.(l) : 

ffj = fqCOS(p q , ft = fqSinipg (7) 

= -fs-simps , ft = f s cos(p s (8) 

As we see that in each basis there are two angles and in [2] it was shown that, to a very 
good accuracy, the mixing can described in terms of single angle f since \ip s — <p q \/\ip s + 
ip q \ « 1 which is also confirmed by a QCD sum rules calculation [5]. 
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In present work we will follow to the first approach, neglecting the mixing angles #o and 
#8 due to their smallness and calculate the coupling constant of p — > 777 decay in framework 
of light cone QCD sum rules (more about of light cone QCD sum rules and its applications 
can be found in [6] and [7]). 

The paper is organized as follows: In section 2 we derive light cone QCD sum rules for 
p — > 777 decay constant, section 3 is devoted to the numerical analysis of the sum rules and 
contain our conclusions. 



2 Light Cone QCD Sum Rules for p — > 777 coupling 
constant 

In this section we calculate the coupling constant of p — > 777 decay using light cone QCD 
sum rules method. In order to calculate this coupling constant we consider the following 
two point correlation function 

U^ = iJ dW** (0\T{j2(x)JftO)} |0) 7 (9) 

where 7 denotes the external electromagnetic field, and JJj and J? are the interpolating 
currents with 7/ and p meson quantum numbers. Here we would make the following remark. 
As we already noted that both of the mixing angles in Eq. (5) are small, in the following 
discussion we will neglect the mixing, i.e we take JJl = Jf . 

At the phenomenological level the Eq.(9) can be expressed as: 

n V <0|^l^(P2)>^(P2)|p(Pl)> 7 (p(Pl)|^|0> nm 

V ^ (Pi ~ rn*)(pl - ml) [W) 

where pi = p 2 + q and q is the photon momentum. The matrix elements entering Eq.(lO) 
are defined as 

<0|J» = m p f p e^ (11) 

(0\4MP2)) = -if*P2» (12) 

The remaining matrix element (i]{P2)\p(Pi))y which describes the Ml transition, can be 
parameterized as, 

(v(P2)\p(pi)) 1 = ee^e^elq/sFiq 2 ) (13) 

where e 1 is the photon polarization vector. Since the photon is real, we need the value of 
F(q 2 ) only at the point q 2 = 0. We can use an alternative parametrization for the ^777 
vertex: 

g 

Lint = 9pn^uual3(d u p u - d u p fI )(d a Ap - OpA a ) (14) 

Tflp 

Comparing Eqns. (13) and (14) we see that 

F(g 2 = 0) = iB!!l (15) 
771 p 
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Using Eqns. (10) - (14), for the physical part of the sum rules we get 

{pl - rrfyipl - m 2 p ) 



Our next task is the calculation of correlator Eq.(l) from the QCD side. The correlator 
receives both perturbative and non-perturbative contributions. In calculation of the non- 
perturbative contributions by the OPE on the light cone one needs to know the matrix 
elements of nonlocal operators between vacuum and the photon states;i.e. ( , y{q)\q^iq\0) 
where Tj is an arbitrary Dirac matrix. These matrix elements can be expressed in terms of 
photon wave functions with definite twist. In calculations we neglect twist 3 three particle 
photon wave functions since their contributions are small. Twist two and twist four photon 
wave functions are defined as [8, 9, 10]: 

(7(9)|9<7a/J5|0) = ie q {qq) f due luqx {(e a qp - epq a )\xw{u) + x 2 [ gi (u) - g 2 (u)}} + 

Jo 

[qx(e a x p - £px a ) + ex(x a qp - q a xp)}g 2 (u)} (17) 
(l(q)\qial 5 q\0) = f -e q e, aPp e a q p x p £ due iU ^(u) (18) 

where for simplicity we use qTq to denote Jf^ = -^{uYu + dYd — 2sYs). 

In Eqs.(17) and (18) e q is the corresponding quark charge, \ is the magnetic suscep- 
tibility, <f(u), and ip(u) are the leading twist two and gi(u), and #2 (if) are the twist four 
photon wave functions. 

After standard calculations we get the following expression for correlator from QCD 
side in the coordinate representation: 

1 x c 

IV = "7i2 16^" ~~ ^ J d4xeW2X {[ Q ( £x ) e c*P^(la x f3 ~ Q(qx)e af3ufl e a Xf3 + 6x 2 e a p up e a qp 
-\2x v e a p ppL e a qpx p \ + fit 2 x 2 [x„e a p pli e a qpx p + 

x^e a p pv e a qpx p \ [ due iuqx ^{u)} /ttV (19) 

J 

The sum rules for g mi are obtained by equating the phenomenological and theoretical 
parts (in Eq.(19) it is necessary to perform Fourier transformation first) of the correlator. 

Performing double Borel transformation on variables p\ = p 2 and p\ = (p + q) 2 on both 
sides of the correlation function in order to suppress the contributions of the continuum 
and higher states (procedure of subtraction in light cone sum rules one can find in [11, 12]), 
and also to remove the subtraction terms in the dispersion relation, we obtain the following 
sum rules for g prjl coupling constant 

g m = 2V3n (1 " Uo){ ^ m2e ^o/m 2 ) + Mu )} (20) 

where Eq(sq/M 2 ) = 1 — e ^ s «/M 2 - 1S ^ e f unc tion usec [ to subtract continuum, s is the 
continuum and 

M 2 2 _ M 2 M 2 

U ° ~ M 2 + Mf ' ~ M 2 + M 2 (21) 
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where M\ and M| are the Borel mass parameters in p and rj channels. Note that in Eq.(21) 
we take into account = Fjj. The masses of r] meson is close to the p meson mass. For 
this reason, it is natural to set Mf = Mf = 2M 2 from which it follows that uo = 1/2. 

3 Numerical Analysis 

In this section we present our numerical result on g m coupling constant. From sum rules 
Eq.(20) we see that for estimating g mi coupling constant first of all one needs to know the 
photon wave function ip(u). It was shown in [9, 10] that the photon wave function do not 
deviate remarkably from its asymptotic form which is given by ip{u) = 1 [9, 8]. The values 
of the other constants appearing in the sum rules are: m p = 0.77 GeV, m v = 0.55 GeV, f = 
0.028GeV 2 . Leptonic decay constant of p meson f p = 0.15 GeV follows from experimental 
result of the p — > e + e~ decay , T(p — > e + e~) = (6, 85 ±0, ll)lfeV[13]. More recent analysis 
shows that the coupling of rj meson with the octet axial vector current is /® = 0.159 GeV 
[2] and this result we will use in our analysis. 

In Fig. 1 we present the dependence of the coupling constant g pm on the Borel parameter 
M 2 at three different values of the continuum threshold: s = 1.4 GeV 2 , 1.6 GeV 2 , 1.8 GeV 2 . 
Since the Borel mass M 2 is an auxiliary parameter and the physical quantities should not 
depend on it, we must look for the region where g pm is practically independent of M 2 . We 
obtain that this condition is satisfied when 1 GeV 2 < M 2 < 1.4 GeV 2 . From this figure 
we also obtained that the variation of s from s — 1-4 GeV 2 to s — 1-8 GeV 2 causes a 
change on the result on g mi of about 10%. Therefore one can say that the result g pm is 
insensitive to s and M 2 . Our final prediction on the coupling constant is 

^ = (1.4 ±0.2) (22) 

where the error is attributed to the variation of s , M 2 and neglected twist three photon 
wave functions. 

At the end we would like to compare our prediction on g mi with experimental result. 
The decay width of the p — > 777 decay is given by 

r(p-^) = ^||lm p (l-mj/mj) 3 (23) 

The experimental value is r(p° — > 777) = (57±10) lfeV[13]. Using this value of r(p° — > 777), 
the g pm coupling constant is obtained from Eq.(23) as: 

9pvi = (1.42 ±0.12) (24) 

which is very close to the sum rule prediction. 

Finally we note that the coupling constant g wrn can be obtained from g prjl with the help 
of the relation g prjl = 3g wrn 
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Figure 1: The dependence of the g m coupling constant on the Borel parameter M 2 at 
three different values of the continuum threshold s = 1-4 GeV 2 , s = 1.6 GeV 2 and 
s = 1.8 GeV 2 . 
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